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Numerical reconstruction of photon-number statistics from photocounting statistics: 

Regularization of an ill-posed problem 
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We demonstrate a practical possibility of loss compensation in measured photocounting statistics 
in the presence of dark counts and background radiation noise. It is shown that satisfactory results 
are obtained even in the case of low detection efficiency and large experimental errors. 
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Photoelectric detection of quantum light [HQ is a ba- 
sic experimental technique in a variety of fundamental 
and applied investigations. In principle, the photon- 
number resolved detectors enable one to determine the 
number of photons in radiation fields. In practice, the 
number of photocounts may significantly differ from the 
number of photons due to losses, dark counts, and back- 
ground radiation. The modern technologies enable one 
to get the detection efficiency near 0.9 and even more Q. 
However, such an improvement leads, as rule, to an in- 
crease in the dark count rate Furthermore, different 
losses occur at all the stages of generation, manipulation, 
and transmission of quantum light. 

The effects of losses and noise in photocounting statis- 
tics can be compensated in two different ways. First, 
the active compensation can be realized by using the ho- 
modyne detection and preamplification of the signal by 
a degenerate parametric amplifier Another possibil- 
ity is a numerical manipulation with the measured data 
- the corresponding technique of loss compensation has 
been discussed in Ref. . The problem is that the corre- 
sponding scries can diverge for small efficiency rj in many 
important cases that require application of a special tech- 
nique of the analytical continuation 0] • 

However, the most serious problem is that the method 
of numerical compensation occurs to be unstable with 
respect to small experimental inaccuracies. Small exper- 
imental errors in photocounting statistics may lead to 
large errors in the reconstructed photon-number statis- 
tics even for large detection efficiencies. This problem 
should be resolved by application of special regularization 
methods. For example, the photon-number statistics of a 
laser radiation has been reconstructed by the method of 
maximum entropy in Ref. The least-squares regular- 
ization technique for loss compensation in photocounting 
statistics has been recently considered in the context of 
the tomography of the quantum detectors 0|. The nu- 
merical compensation of losses in multi-pixel detectors 
has been discussed in llOj. The method of maximum- 



likelihood estimation demonstrates satisfactory results 
for quantum-state reconstruction in the presence of losses 
ll| and, consequently, it can be applied for loss compen- 
sation in photocounting statistics. An alternative tech- 
nique of the regularization, which requires measurements 
with different values of the efficiency, has been proposed 
in Ref. \v^. 

In the present contribution we reexamine the method 
of numerical compensation of losses. We demonstrate 
that the regularization of the corresponding ill-posed 
problem (see, e.g., [3] and [3]) enables one to apply 
this technique even for low values of the efficiency r\. 
Moreover, our consideration includes numerical compen- 
sation of dark counts and effects of background radia- 
tion. Besides, the method demonstrates satisfactory re- 
sults under a realistic assumption that the efficiency and 
the noise-counts rate are known with a certain inaccu- 
racy. Wc demonstrate that the proposed technique can 
be applied for different photocounting statistics including 
highly-nonclassical cases. 

Let us start with consideration of a single-mode quan- 
tum light characterized by the density operator g. If n 
is the corresponding photon- number operator and \n) is 
its eigenvector, the photon-number distribution is given 
by " 



cigen 



(n\Q\n) 



Tr 



— exp (-n) : g 



(1) 



In the presence of losses, dark counts, and background 
radiation it differs from the photocounting distribu- 

urn. 



tion 



T n = Tr 



(77 h +N nc ) 7 



exp (-7771 - 2V nc ) : g , (2) 



where rj and A^ nc are the efficiency and the mean num- 
ber of noise counts, respectively. The aim of this work 
is to develop a mathematical technique for reconstruc- 
tion of the photon-number distribution p n from the 
experimentally-measured photocounting distribution V n - 
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The photocounting distribution V n is expressed in 
terms of the photon-number distribution p n as (see [E 

G3) 



Pro = X] ^"'l™ ("> ^> c ) 
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Here, 



S m \ n (V, N nc ) = e- N -NZ- n V n ~, L™- f^hJl 

m! \ V 



(4) 



for m > n and 



S m |n (n, A nc ) = e -"«(l-7 ? )»-"y i» 



A nc (n - 1) 
V 



(5) 

for m < n are the probabilities to get m photocounts 
under the condition that n photons are present. L™(x) 
is the Lagucrrc polynomial. Expression ((3J is formally 
inverted as 
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where 

S n\L 0?> N nc ) = (n + 1, n - m + 1; 
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S-^ („, ]V nc ) = e N ™$ L + l, m - „ + 1; N ^ 
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for m > n is the matrix inverse to S m \ n (rj, N nc ). 
$ (n, m; a;) is the Kummcr hypcrgeometric function. 

As mentioned, expression j6j as the solution of Eq. (|3|) 
is unstable with respect to small experimental inaccu- 
racies of P n . Moreover, similar to the case of zero noise 
counts [E0] i this series diverges in many important cases. 
Hence, Eqs. ([6|)-([8|) cannot be applied in the general case. 
This fact is a consequence of a mo re g eneral statement 
that Eq. j3]) is an ill-posed problem [13| . Such a problem 
can be treated by using the appropriated regularization 
methods [IE Hull. 

A priori information, such as numbers at which 
photon-number and photocounting distributions can be 
truncated [la], is used in the regularization of the ill- 
posed problem. In addition, the basic properties of the 
photon-number distribution, 



Pn > 0, 
5>n = 1, 



(9) 
(10) 



arc applied in the considered case. Other a priori infor- 
mation can also be useful for the regularization of the 
ill-posed problem depending on a given physical situa- 
tion. 

In different problems of quantum optics the least- 
squares inversion and the Tikhonov regularization lead 
to satisfactory results (for a review see, e.g., [Hj] ). In 
this contribution, we apply the Landweber algorithm |l9l| 
adopted to the regularization of similar problems [17j - 
a technique, which demonstrates a good computer com- 
patibility [2Q]. The projected Landweber algorithm [2lj| 
is the iteration process, 



Si) 



+ x(s^v-s^s P ^- 1 ^] , (11) 



where pV) = ^p\P j is the jth iteration for the photon- 
number distribution, V = {V m }, S = {S m \ n }, and x 
is the relaxation parameter. 11c is the projector on the 
closed convex set C defined by Eq. j9]) and, in special 
cases, by other additional conditions. Condition (|10[) can 
be used to track the accuracy of the obtained results. The 
starting values are usually chosen as p„^ = 0. 

To illustrate the method let us give some numerical 
simulations. We start from the thermal state, 



1 



Q 



1 + n t h V 1 + nth 



(12) 



with n t h = 30, and derive from Eq. ^ the photocount- 
ing distribution V n for n = 0.34 and N nc = 0.30. The 
measured data are simulated with i/ = 5x 10 4 sampling 
events. The corresponding relative error (in terms of 
the Euclidian norm) is Sp = 0.03. The simulated data 
are then used as an input of Landweber algorithm (lll|) 
for the reconstruction of the photon-number distribution 
with inaccurate values of the efficiency f) = 0.35 and the 
mean number of noise counts N nc — 0.29. The result of 
this procedure (see Fig. [1]) is in a reasonable agreement 
with the initially chosen photon-number distribution, the 



relative error is S n 



0.05, and the relative residual is 



dp = 0.019. It is worth noting that in the given example 
seriesj3|) diverges even in the absence of experimental er- 
rors [a]. Nevertheless, algorithm (fTTj) demonstrates high 
efficiency for the considered case. 

Another example is the single-photon-added thermal 
state (SPATS), 



rath (1 + nth) \l+rath 



rath 



(13) 



with nth = 10- Such a state has been recently realized 
experimentally and its nonclassical properties have been 
verified [22| . The numerical simulation is performed for 
v = 5 x 10 3 , n = 0.7764, and A nc = 0.748, with the rel- 
ative error Sp = 0.025. The photon-number distribution 
(reconstructed with fj = 0.77, A nc = 0.75) is shown in 
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FIG. 1: (Color online) The photon-number, p n , and photo- 
counting, V n , distributions of the thermal state, nth = 30. 
The circles and triangles are the initially chosen and recon- 
structed (fj = 0.35, Nxxc = 0.29) photon-number distributions, 
respectively. The asterisks show the simulated photocounting 
distribution (y = 5 x 10 4 , rj = 0.34, N nc = 0.30). 




4 6 8 

Number of photons 

FIG. 3: (Color online) The photon-number distribution, p n , 
of the SPATS, nth = 10. The circles and triangles are the 
initially chosen and reconstructed (y — 5 x 10 5 , rj — 0.7764, 
N nc = 0.748) distributions in Eq. ([BJ, respectively. The result 
of such a reconstruction demonstrates much stronger noise ef- 
fect in comparison with the Landweber algorithm (see Fig. [2J . 



Fig. [21 The relative error is b p = 0.041 and the relative 
residual is b p = 0.020. In this example series ([6]) formally 
converges. However, attempts to apply this series for the 
direct reconstruction of p n result in the large noise effects 
caused by small experimental inaccuracies even for suf- 
ficiently large number of sampling events (y = 5 x 10 5 , 
bp = 0.0079) and exact values of rj and N nc (see Fig. [3]). 
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FIG. 2: (Color online) The photon-number, p n , and photo- 
counting, Vn, distributions of the SPATS, nth = 10. The cir- 
cles and triangles are the initially chosen and reconstructed 
(f) — 0.77, Nnc = 0.75) photon-number distributions, respec- 
tively. The asterisks show the simulated photocounting dis- 
tribution (y = 5 x 10 3 , f) = 0.7764, iV nc = 0.748). 



Algorithm (fTTj) demonstrates a reasonable agreement 
with the initially chosen photocounting distribution even 
in the case of rather large error in the simulated data. 



To illustrate this fact, we consider the superposition of 
the coherent states, 



1 



2(l + e- 2 M 2 



= (|a) + |-a>), (14) 



\a\ = 23.9. The initial data are simulated with v = 
5 x 10 3 , <q = 0.613749, and N nc = 1.763442. The relative 
error is bp = 0.098. The photon-number distribution 
is reconstructed for r) = 0.59, N nc = 1.77 (see Fig. 
The relative error of the reconstructed distribution, b p = 
0.125, and the relative residual, b p — 0.051, are of the 
same order as for the initially simulated data. It should 
be stressed that for the state (fT4")) p n ^ only for even 
photon numbers n. This a priori information is used in 
the projector lie, [cf. Eq. (fTTj) ]. 

In conclusion, we have obtained the expression for pho- 
tocounting distribution in terms of the photon-number 
distribution. However, the inverted expression cannot 
be used in the most practical situations - it is unstable 
with respect to small experimental inaccuracies and the 
corresponding scries can diverge. At the same time, the 
regularization of this ill-posed problem by the Landwe- 
ber algorithm enables one to compensate losses and noise 
counts even for low efficiencies, high noise-counts rates, 
and inaccurate knowledge of their values. 
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FIG. 4: (Color online) The photon-number, p n , and photo- 
counting, Vn, distributions of the superposition of the co- 
herent states |q| 2 = 23.9. The circles and triangles are the 
initially chosen and reconstructed ( fj — 0.59, N nc = 1.77) 
photon-number distributions, respectively. The asterisks 
show the simulated photocounting distribution (y — 5 x 10 3 , 
r, = 0.613749, iVnc = 1.763442). 
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